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ABSTRACT 

The baryon acoustic oscillations are a promising route to the precision measure of the cosmological 
distance scale and hence the measurement of the time evolution of dark energy. We show that the non- 
linear degradation of the acoustic signature in the correlations of low-redshift galaxies is a correctable 
process. By suitable reconstruction of the linear density field, one can sharpen the acoustic peak in 
the correlation function or, equivalently, restore the higher harmonics of the oscillations in the power 
spectrum. With this, one can achieve better measurements of the acoustic scale for a given survey 
volume. Reconstruction is particularly effective at low redshift, where the non-linearities are worse 
but where the dark energy density is highest. At z = 0.3, we find that one can reduce the sample 
variance error bar on the acoustic scale by at least a factor of 2 and in principle by nearly a factor of 
4. We discuss the significant implications our results have for the design of galaxy surveys aimed at 
measuring the distance scale through the acoustic peak. 

Subject headings: large-scale structure of the universe — distance scale — cosmological parameters — 
cosmic microwave background 
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1. INTRODUCTION 

The late-ti me acceleration o f the expansion rate of 
the Universe ('Riess et al.l ll998t iPerlmutter et aP 11999 ^: 
^css ct al. 2004) argues for a remarkable change to 
our understanding of the forces of nature. Choos- 
ing betwen the exotic explanations to this surpris- 
ing phenomenon may be possible through precision 
measure of the expansion rate of the Universe and 
of the growth of cosmological structure over time. 
Sound waves propagating in the first 400,000 years 
after the Big Bang produce a characteristic length 
scale in the anisotropics of the mi crowave background 



and in the clustering of galaxies ('Peebles 
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IIoltzmann ^l989": Hu & Sugivama 1996; Eiscn stein fc Hul 
19981) . CMB experiments have clearly detected these 
fluctuations in both the temperature and polarization 
power spectru m (for a summary of recent results, see 
iSpergel et al.l l2006). With large galaxy surveys, we can 
detect this acoustic signature as a peak in the cor- 
relation function at ^ 150 Mpc or as a ha rmonic se- 
quen c e of oscillations in the power spectrum IjCole et al.l 
120051 lEisenstein et al.ll2005D . This length scale can be 
measured as a characteristic angle on the sky, yield- 
ing the angular diameter distance as a function of red- 
shift, and as a characteristic difference in redshift of 
galaxy pairs alon g the line of sight , yielding the Hubble 
parameter H( z) l|E isensteinl l2003t iBlake fc Glazebrookl 
|2003; Hu & H aimanI 120031 iSeo fc EisensteinI 12003(1 . Be- 
cause of the very large scale, the acoustic signature re- 
mains in the linear regime even today and is therefore a 
highly robus t method for measur i ng the cosmological dis- 
tance scale llMeiksin et al.lll999t iSeo fc Eisensteinl 120011 
lEisenstein et a,]Jl200(i(l. 

^ Steward Observatory, University of Arizona, 933 N. Cherry 
Ave., Tucson, AZ 85121 

^ Dept. of Astrophysical Sciences, Peyton Hall, Princeton Uni- 
versity, Princeton, NJ 08544-1001 



However, as the universe evolves, the acoustic signa- 
ture in the correlation function is broadened by non- 
linear effects. Non-linear gravitational structure forma- 
tion and redshift distortions move galaxies away from 
their original locations, blurring out the peak at 150 Mpc 
separation, or, equivalently, erasing the higher harmon- 
ics in the power spectr um that represent smaller scales 
IjEisenstein et aljf2fl0f]tl . Even though the baryon bump 
remains measurable in the correlation function, these 
nonlinear effects reduce the precision of the distance scale 
that can be measured from a given volume of space, by 
roughly a factor of three at the present day. This loss is 
significant because survey volume is precious: larger sur- 
veys are more expensive and the amount of volume that 
we can survey in the low-redshift Universe is limited. 

We argue in this paper that this loss of precision is 
avoidable. The blurring of the acoustic peak is largely 
due to bulk flows and super-cluster formation, effects 
that are generated by gravitational forces on large scales. 
The same map of galaxies intended to measure the acous- 
tic scale is an accurate map of the large-scale gravita- 
tional source terms. One can essentially run the gravita- 
tional flow backwards to restore the acoustic peak nearly 
to its linear regime shape. 

2. NON-LINEARITIES 

Numerical simulations show that the advancing scale 
of non-linear gravitational collapse e rases the higher har- 
monics of the acoustic oscillations llMeik sin et al ]ll999l: 






geo fc Eisnnstnin'2005':'Snrin ge] et a,1J2005HWhitfl2005r 
Eiscnstein et al. (2005, 2006) note that the harmonic se- 
quence in the power spectrum correspond to a single 
peak in the correlation function and that the damping 
envelope corresponds to the broadening of this peak. In 
the configuration-space view, the source of the late-time 
broadening is clear: matter is b eing moved by r oughly 
10 Mpc from its initial position. lEisenstein et al.l ((200(tl 
build a model for the non-linearity in terms of the dif- 
ferential motion of pairs initially separated by 150 Mpc. 
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Fig. 1. — The cumulative variance in the diflFerential motion 
of pairs initially separated by 150 Mpc as a function of cutoff 
wavenumber, normahzed to the total variance. These curves are 
calculated in the Zel'dovich approximation according to the formu- 
lae in^isenstein et al. (2006). The displacement along the sepa- 
ration vector is shown as the solid line; the displacement along a 
single direction perpendicular to the separation vector is shown as 
the dashed line. Most of the integral is supported by wavenunibers 
between 0.02 and 0.2/iMpc~^. The horizontal dashed line is drawn 
at 75%, where half of the rms displacement has been fixed. 



The final large-scale correlation function is simply the 
convolution of the linear correlation function with the 
distribution of differential motions. 

One can ask what scales are responsible for the dif- 
ferential motion . Using the Zel'dovich approximation 
l)Zeldovichlll970(l . the second moment of the distribution 
of the differential motion can be written as an integral 
over the power spectrum (Eisenstein ct al. 2006). The 
cumulant of the integrand is plotted in Figure ^ Here 
we see that most of the motion is generated at A: « 
0.1/iMpc~^. Wavenumbers smaller than 0.02/iMpc~^ 
contribute little, because these perturbations affect both 
points equally. Wavenumbers larger than 0.2/iMpc~^ 
contribute less because the CDM power spectrum is fairly 
red on these scales. 

The large-scale velocity field is responsible for most of 
the nonlinear effects that appear to "erase" the high or- 
der acoustic peaks in the power spectrum. The motions 
are dominantly due to bulk flows and the formation of 
superclusters. Small-scale formation of halos is subdomi- 
nant; even big halos only pull material from an average of 
h(Mlm^^hr^ M0)i/3/i-i Mpc away, about half the bulk 
flow motions. 

3. RECONSTRUCTION 

Our major point is that because the scales of inter- 
est are large, the motions of galaxies can be well mod- 
eled by perturbation theory and in principle can be mea- 
sured and removed. In the simplest terms, the bulk flows 
are generated by exactly the density perturbations that 
are being surveyed to measure the acoustic oscillations. 
The connection between the density and velocity field on 
these large scales is nearly that of linear theory, so one 
can predict the velocity field and undo the motion of the 
galaxies. 

In more detail, reconstruction of the velocity field 
or the linear density field is a su bject with consider- 
able history. Peebles (,19891 Il990() pioneered the sub- 
ject by reconstructing the trajectories of Local Group 



galaxies using the principle of least action, with the 
goal of constraining the local val ue of J?n- The sim- 
ple "Gaussianization" method of iWeinbcrg* il992() ap- 
proac hed the problem from a cosmological perspec- 
tive. "Nusscr fc Dekell |1992) and G ramannI (1993) re- 
constructed primordial density fields with a technique 
based on the success of the Zeldovich approximation in 
the quasilinear regime, similar to our technique described 
below. Naravanan & Croft (1999) developed these tech- 
niques further (and offered a comparison of several differ- 
ent methods). Monaco & Efstathiou (1999) investigated 
a self-consistent, but more complicated, iterative scheme 
in which the initial densities source the Lagrangian map, 
rather than the final densities. Pursuit of this recon- 
struction method further may warrant using higher or- 
der relati onships be tween pe culiar velocities and densi- 
ties (e.g.. iChodorow ski et al] ll998tl. The technique s of 
iCroft fc Gaztanaeal (,1997.) and iBrenier et all 1)2003(1 are 
based on the principle of least action again, with the ap- 
proximation that particles move on straight trajectories. 
As far as we know no one has attempted a full cosmo- 
logically relevant proper least action integration back- 
wards in time, as the computational resources would be 
prohibitively expensive. Reversing the sign of gravity 
and running an evolved simulation backwards would not 
work, at least in the realistic case in which there are er- 
rors or unconstrained dimensions in the galaxies' phase 
space, because the decaying modes would overpower the 
growing modes as one went back in time. 

Restoring in full the acoustic signature at fc < 
0.2/1 Mpc^^ is an undemanding application of these re- 
construction techniques. To demonstrate that recon- 
struction can help, we here show a simple method. 
We take the 2 = 0.3 output s from the 512 hr^ Mpc 
box N-body simulations from iSeo fc Eisensteiia l(200,'Tf l 
{N = 256"^ particles), compute the density field, Fourier 
transform and filter with a Gaussian of lQh~^ Mpc or 
20/1^^ Mpc width. From this, we predict the linear the- 
ory motion using V ■ q = —6, where qis the Lagrangian 
displacement field and S is the fractional overdensity. 
We then move the particles by —q. We do the same 
for a reference grid of smoothly distributed particles. In 
redshi ft space, we account for the linear redshift distor- 
tions (|Kaisei<ll987j) by boosting the displacements along 
the line of sight by an additional factor of 1 -I- /, where 
/ = d{\n D) / d{\n a) , D is the linear growth function, and 
a is the scale factor of the Universe, A new density field 
is defined by the difference of the density field of the real 
particles and that of the re ference particles. Note that, 
in contrast to the work of iSirko fc Soergeil (|2006D . the 
point is not to move all of the particles back to their ini- 
tial location, but rather to move the measured densities 
back to their initial location. 

The power spectrum of this field is shown in Figure 
121 This Figure is the average of power spectra from 30 
simulations at z — 0.3. One can see the degradation 
of the higher harmonics at z = 0.3 in the uncorrected 
density field, compared to the initial spectrum at z — 
49, and see that the reconstruction has partially restored 
them. Figures El and ^ show the correlation functions 
in real and redshift space. The acoustic peak is nearly 
fully restored in real space. In redshift space, the peak 
is considerably improved but not fully fixed. 
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Fig. 2. — The matter power spectrum after reconstruction by 
the hnear-theory density-velocity relation, with the density field 
Gaussian filtered. The bottom panel shows the real-space power 
spectrum; the top panel shows the spherically averaged rcdshift- 
space power spectrum. The black solid line shows the input power 
spectrum at 2 = 49; it has been displaced in the top panel for 
clarity. The blue short-dashed line shows the matter power spec- 
trum at 2: = 0.3; one can see that acoustic peaks have been lost. 
In the bottom panel, the red dot-dashed line and magenta long- 
dashed line show the effects of reconstruction for 20h~^ Mpc and 
10h~^ Mpc Gaussian filtering, respectively. In the top panel, both 
lines show 10h~^ Mpc filtering; the red dot-dashed line is without 
finger of God compression, while the magenta long-dashed line in- 
cludes compression. The increase of power at large wavenumbers 
is essentially irrelevant to the quality of the acoustic signature; one 
would in practice marginalize over these broadband changes. 



When one predicts the large-scale displacement field, 
one is also predicting the large-scale velocity field 
and hence the correction for large-scale peculiar ve- 
locities in reds hift space. This is important because 
lEisenstein et al . (2006) find that redshift space distor- 
tions degrade the radial measurement of the acoustic 
peak. On large scales, the real-space displacements of 
particles are in the same direction as their peculiar ve- 
locity distortion, so the degradation of the acoustic peak 
is worse in redshift space than in real space. Reconstruc- 
tion can fix this. 

However, there are also redshift-space distortions from 
small-scale peculiar velocities, i.e., fingers of God. Clus- 
ters of galaxies appear as long cigars along the line of 
sight in redshift space. For the purposes of determining 
bulk flows, one should simply compress these fingers of 
God back to some approximation of their real-space loca- 
tion. Without this step, the fingers of God get stretched 
out further by the reconstruction, degrading the acous- 
tic signature. To show that finger-of-God compression 
can help the reconstruction, we have identified clusters 
in redshift space with an anisotropic friends-of-friends al- 
gorithm and moved all cluster particles to the center of 
mass of the cluster. The correlation function that results 
from running our reconstruction on the compressed den- 
sity field is shown in Figure ^ One sees a modest but 
useful improvement. 

The simple reconstruction described above has not 
fully restored the linear acoustic scale, particularly when 
beginning from redshift space. We expect that more so- 
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Fig. 3. — The real-space matter correlation function after recon- 
struction by the linear-theory density-velocity relation, with the 
density field Gaussian filtered. The black solid line shows the cor- 
relation function at 2 = 49. The blue short-dashed line shows it 
at 2 = 0.3; the acoustic peak has been smeared out. The red 
dot-dashed and magenta long-dashed lines show the effects of re- 
construction for 20h~^ Mpc and 10h~^ Mpc Gaussian filtering, re- 
spectively. Even this very simple reconstruction recovers nearly all 
of the linear acoustic peak. 
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Fig. 4. — The redshift-space matter correlation function after 
reconstruction by the linear-theory density-velocity relation, with 
the density field Gaussian filtered. The black solid line shows the 
correlation function at 2 = 49. The blue short-dashed line shows 
the redshift-space correlation function at 2 = 0.3; the acoustic peak 
has been smeared out. The black dotted line shows the real-space 
correlation function for comparison. The red dot-dashed line line 
shows the effects of reconstruction for a 10h~^ Mpc Gaussian filter- 
ing; the magenta long-dashed line is the result when one compresses 
the fingers of God prior to the reconstruction. These reconstruc- 
tions significantly improve the acoustic peak. 



phisticated reconstruction methods will produce further 
improvements. The small end of the range of the scales of 
interest are in the quasi-linear regime, and our assump- 
tion of linear theory for both the continuity equation and 
the redshift distortions is only a first approximation here. 
We use a Fisher matrix calculation to estimate how 
much the reconstruction has improved the recovery of 
the ac oustic scale. Our calcula tion is based on the meth- 
ods in lSeo fc EisenstemI l)2003j) but with the derivatives 
multiplied by a Gaussian filtering that i s tuned to match 
the p air-wise Lagrangian displacement IjEisenstein et alj 
|2006|) and that visually reproduces the smearing of the 
acoustic peak. We focus here on the spherically-averaged 
acoustic scale; we will present anisotropic results (i.e., 
separate estimates for the angular diameter distance and 
the Hubble parameter) in a future paper. We find that 



at z — 0.3, in the absence of reconstruction, a survey 
of Ih^^ Gpc'^ should produce a distance measurement of 
1.4% from real space and 1.9% from redshift space. The 
real-space value is in good agreement from the jack-knife 
estimates from the simulated real-space power spectra in 
(jSeo fc Eisenstein 2005') . The simple reconstruction pre- 
sented in this section improves these measurements to 
0.75% and 0.95% in real and redshift space, respectively. 
This is a factor of two improvement. We find that if the 
reconstruction were perfect, one could achieve 0.5% dis- 
tances (although shot noise from reasonable galaxy sam- 
ples would degrade this to 0.55-0.60%). Hence, the raw 
precision at z = 0.3 is roughly a factor of 3.5 worse than 
the true cosmic variance, and the simplest reconstruction 
returns the first factor of two. At higher redshifts, the 
cosmic variance per unit volume won't change, but the 
data sets with or without reconstruction will be closer to 
that ideal. 

4. DISCUSSION 

We have demonstrated that density-field reconstruc- 
tion can restore the linear-regime contrast of the baryon 
acoustic signature. Physically, the galaxy map locates 
the superclusters and voids and permits one to predict 
the large-scale flows resulting from these objects. By 
moving objects back to their initial location, one can cor- 
rect the effect that these flows have on the characteristic 
separation of galaxies produced by the acoustic waves in 
the early universe. By restoring the linear-regime cluster- 
ing, one can improve the precision of the measurement of 
the acoustic scale from galaxy redshift surveys and hence 
the available constraints on the cosmic distance scale and 
dark energy. 

The opportunity for density reconstruction alters the 
optimization of surveys for baryon acoustic oscillations. 
Low-redshift surveys were thought to suffer in perfor- 
mance per unit volume; with reconstruction, this is not 
the case. Given that dark energy is more important (and 
hence easier to measure) at lower redshift, this will tend 
to push the optimal redshift range lower. On the other 
hand, lower redshift surveys were generally designed at 
lower target densities, because the suppression of the 
higher harmonics meant that it was not necessary to mea- 
sure k « 0.2/iMpc~^ well. Now the opportunity exists to 
get acoustic oscillation information from the higher har- 
monics, so that one might want to reduce shot noise and 
measure k « 0.2h Mpc~ well. In practice, with multi- 
object spectrographs, the number density at the highest 
redshifts will drop due to the flux limits imposed by a 
given exposure time, but at lower redshifts one can use 
the extra galaxies available to that flux limit, if one has 
sufficient number of fibers or slits. 

Do the redshift surveys aimed at acoustic oscillations 
have enough information to do density-field reconstruc- 
tion? Acoustic oscillation surveys are designed to bal- 
ance shot noise and sample variance at the wavenum- 
bers where the acoustic oscillations are found, namely 
0.1 — 0.2/iMpc~^. If a survey is measuring at fc < 
0.2/iMpc~^ with shot noise below sample variance, then 
it is producing a reasonable fidelity map at exactly the 
wavcnumbcrs required to do the reconstruction, as shown 
in Figure ^ Turning this around, if the map is insuf- 
ficient for reconstruction, it will also be too noisy to 



measure the power spectrum well (i.e., approaching the 
sample variance limit) at the higher harmonics that one 
was hoping to improve. Indeed, within the approxima- 
tion of the linear-theory density-displacement relation, 
the desired weighting of each Fourier mode to be used 
in reconstruction is nP/{l -t- nP), where n is the co- 
moving number density and P is the power spectrum 
at that wavenumber. This is the same familiar factor 
that deterr nines the measurement of modes in the power 
spect rum l|Feldman . Kaiser, fc Peacock! 119941 : iTeemarfl 
Il997|) . Hence, surveys designed to measure the acous- 
tic oscillations with moderate signal-to-noise ratio maps 
will be well suited to reconstruction. 

One must remember that the reconstruction need not 
be perfect. The propagation of sound waves at the epoch 
of recombination has mild dispersion that gives the linear 
regime acoustic peak a width of about 30 Mpc FWHM or 
8/i~^ Mpc rms width, comparable to the pairwise displac- 
ment s of galaxies even at low redshifts l|Eisenstein et alj 
I200(J) . Once we reduce the errors below this characteris- 
tic width, there is little improvement in the acoustic os- 
cillations measurements. Thus, reconstruction is easier 
at higher redshifts: the raw displacements are smaller, so 
one can accomplish sufficient reconstruction using only 
larger scales. Even at z = 3, it is likely that reconstruc- 
tion will benefit the recovery of the line-o f-sight a c oustic 
peak (and hence H{z)), because Eisenste in et al.l l)200(Tt) 
shows that the redshift-spacc displacements are nearly 
double the real-space ones at z > 1 and hence are not 
quite negligible even at z = 3. 

Reconstruction places a premium on surveys with con- 
tiguous area. Surveys with lots of gaps on 10-100 Mpc 
scales will not measure the density field well enough to 
do these non-linear corrections: the source of bulk flows 
could be hidden in the gaps. Some gaps, e.g., due to 
bright stars, are inevitable and will degrade the recon- 
struction. Holes of 1 Mpc and smaller are less important; 
these are not biasing the density field at fc = 0.1/iMpc~^ 
much. Of course, the reconstruction techniques them- 
selves must be able to deal with this mildly gappy data. 
Gaps are less of an issue at z > 2 where reconstruction 
is not a large advantage. 

Photometric redshift surveys for acoustic oscillations 
will not measure the 3-dimensional density field on the 
scales required to do density-field reconstruction. This 
is not a problem at z > 2, where the real-space density 
field is not degraded by non-linearities much anyways, 
but is a disadvantage relative to spectroscopic surveys at 
z< 2. 

Much of the work on reconstruction has been focused 
on constructin g the velocity field from the density field, 
or vice versa IjNusser et al.lll99lj) . We are actually in- 
terested in constructing the displacement field from the 
density field. In linear theory, the velocity and displace- 
ment fields are the same, but differences enter at higher 
order in perturbation theory. We expect that much of 
the work on velocity-field reconstruction will be applica- 
ble to the displacement-field problem, and we note that 
in detail the displacement field may be somewhat easier 
to infer, as virialization affects the velocity field substan- 
tially but the displacement field very little. 

The galaxies used in large-volume baryon acoustic os- 
cillation surveys are typically biased, and it will be im- 



portant to build reconstruction techniques that can han- 
dle this bias. Again, one must remember that only mod- 
est performance is required of the reconstruction; if one 
can fix half of the displacements, one will have achieved 
most of the gains. All that is required is that light trace 
mass reasonably well on large scales. Th is is well estab- 
lished by the simple behavior of bias (e.g. lTeemark et aO 
I2004H. the detection of l arge-scale redshift distortions 
(e.g. IPeacock et al.l l200lj) . and the gala xy-mass cross- 
correlation measured from weak lensing IjSheldon et al.l 
12003') . The errors that enter from imperfect knowledge 
of the underlying cosmology or exact level of bias (e.g., 
the value of the redshift distortion factor /9) are much 
smaller than the required reconstruction precision. 

Finally, while reconstruction can improve the preci- 
sion of the acoustic scale measurement, one must also 
verify that the measurement remains unbiased. This 
will require further work, as one must also decide in de- 
tail how one is measuring the distance information from 



the observed clustering, particularly in the case of an 
anisotropic clustering analysis. One would expect recon- 
struction to at least partially fix any bias that enters from 
non-linear gravity, but this has not yet been studied at 
the sub-percent level. 

In conclusion, we believe that density-field recon- 
struction offers a significant opportunity to improve the 
baryon acoustic peak method by a factor of 2 to 3, par- 
ticularly at lower redshifts. 
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